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Abstract 

A criterion to locate tricritical points in phase diagrams is proposed. 
The criterion is formulated in the framework of the Elementary Catas- 
trophe Theory and encompasses all the existing criteria in that it applies 
to systems described by a generally non symmetric free energy which can 
depend on one or more order parameters. We show that a tricritical point 
is given whenever the free energy is not 4-determined. An application to 
smectic-C liquid crystals is briefly discussed. 

1 Introduction 

Despite a common belief, Catastrophe Theory (CT) can provide not only qual- 
itative insight but also quantitative results. This is particularly true for phase 
transitions of systems whose free energy depends on several order parameters, 
such as those occurring in liquid crystals. Usually, CT is applied to concrete 
examples by identifying the elementary catastrophe associated with each transi- 
tion. Calculations can then most easily and efficiently be done on this equivalent 
form. In doing this translation from a physical to a pure mathematical realm, 
one retains the intuition of the problem, but at the same time looses the ability 
to write explicit expressions for the physical quantities involved in the transi- 
tion. 

One of the aims of the present paper is to show how quantitative equations can 
be drawn from CT together with a geometric and intuitive description of their 
meanings when the physical conditions are recast in a diffeomorphic invariant 
form. 

The particular case of tricritical points in 4-order parameter systems is stud- 
ied. Tricritical points are points in thermodynamic phase space where a phase 
transition of a complex kind takes place involving the meeting of a second-order 
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transitions with a line of first-order transitions. In mathematical terms, they 
turn out to be the points with the highest degeneracy in a butterfly Af catastro- 
phe. We then apply our results to a specific example drawn from liquid crystals 
theory, namely smectic-C liquid crystals. A possible subsequent application is 
the study of tricritical points in biaxial liquid crystals. 

This paper is organized as follows. In Section[2]we introduce the mathemat- 
ical background necessary to understand the rest of the paper. In particular 
we omit a few technical details that are important but can be easily studied 
form the general literature and could otherwise obscure the reading of the pa- 
per. Section [3] is devoted to develop the main result: we derive a criterion to 
identify tricritical points in the control space. In Section [4] we specify our main 
result to simpler cases which are physically important and show that our cri- 
terion reduces to the ones already known in the literature when additional ad 
hoc assumptions are made. An application to smectic liquid crystals is given in 
Section [5] and the conclusions are drawn in Section [HI 

2 Mathematical background 

In the present Section we recall the very basic features of CT that will be useful 
in the following. We refer the reader to the excellent introductions [H [51 [5] as 
well as to the review paper [7j for further details. 

Let y(x. A) be a real-valued smooth function with variables x — {xi, X2, . ■ ■ , Xn) G 
M" and parameters A = (Ai, A2, . . . , Xp) G M.P. For each fixed A one obtains a 
function Vx ■ R" — >• R and we are thus enabled to interpret F as a A-parameter 
family of smooth functions from R" to M. The function Vx has a singular point 
at xq for A = Ao if 9x^6^0 (^o) = Vx. We assume that xq = 0, Aq = and 
V(0, 0) = 0, as the general case is deduced from this by translation. We aim at 
analysing the behaviour of Vx in the neighbourhood of 0, as the parameters A 
are varied. Since we are interested in a local study, we lump together all func- 
tions which coincide with Vx near and call this set of functions the germ of Vx 
at 0. For later reference, let m„ denotes the set of germs /(x) with /(O) = 0. 
More generally, we denote by the set of germs in m„ such that all their 
partial derivatives of order less than k vanish at 0. These powers of m„ form 
a descending chain, i.e. m„ D 3 D . . .. It can be shown that m^j is an 
ideal generated by all monomials of homogeneous degree k. 

The germ Vx is said to have a non-degenerate or Morse singular point at 
if its Hessian is there nonsingular, i.e. the rank of the Hessian is n. Otherwise 
the singular point is degenerate or non-Morse. At a Morse singular point, the 
germ Vx exhibits a minimum, a maximum or a saddle point and is structurally 
stable, that is all the germ obtained by a small change of A have the same kind 
of singularity. Moreover, a theorem of Morse guarantees that there exists a 
coordinate transformation which permits to write ^(x, A), in a neighbourhood 
of 0, as 
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(1) Ki(x) ^-{xl + ... + xl) + {xl^, + --- + xl) 

where k is the index of Hessian matrix at 0, i.e. the number of negative eigen- 
values of the matrix. 

At a non-Morse singular point the germ is structurally unstable and varying 
the control parameters it undergoes a catastrophe, that is neighbouring germs 
may possess a different number of singular points with different nature. The 
fundamental example is V\{x) = x^ — Ax, which has one minimum when A > 
and none when A < 0. Non-Morse singularities can be further classified by 
means of corank, determinacy and codimension of Vx at 0. The corank, c, is 
the number of eigenvalues of the Hessian matrix of Vx which are nought. We 
may equivalently think of it as the number of directions in which the germ is 
degenerate. It can be proved that, in a suitable coordinate system, ^(x. A) can 
be written as the sum of a Morse, Vm, and non- Morse part, T^nm, that is 

(2) y(x,A) = l/M(x,) + KM(xe,A) 

where Xe = {xi, . . . , Xc} and = {xc+i, ■ ■ ■ ,Xn} are named essential and 
inessential variables, respectively. This result is called splitting lemma because 
it allows to split the variables into two classes. Unlike the inessential variables, 
the essential variables are involved in the structural instability and the kind of 
catastrophes which can occur depends only on their number, i.e. the corank c 
of the singularity. If the number of control parameters is not greater than five, 
Thom's theorem states that the non-Morse part can be put into a canon- 
ical form, named elementary catastrophe, which is given by the sum of a germ 
which only depends on Xe, said catastrophe germ, and a family of germs which 
depend on Xe as well as on A, said catastrophe perturbation. The latter is also 
named unfolding. The idea is that the degenerate critical point of Vx can be 
potentially unfolded under a perturbation into several non-degenerate critical 
points, which then appear in the nearby germs. Table [1] gives the list of the 
corank one elementary catastrophes. 

Determinacy, k, and codimension, r, of Vx at permit to determine both the 
catastrophe germ and germ to which T^^m is equivalent in a suitable coordinate 
system. We define the k — jet of the germ /(x), which we write as j'^/(x), to 
be the formal Taylor expansion up to including terms of order k. A germ /(x) 
is k-determined if every germ ^(x) having f = j^g is diffeomorphic to /. The 
determinacy of a one- variable germ is trivial since it is simply given by the first 
non-zero term in its Taylor expansion. In more than one dimension this is no 
longer the case. The germ x\x2, for instance, is four-determined yet its Taylor 
expansion contains only three-order terms. It is worth to notice that a germ at 
a Morse singular point is two-determined. 

Algebraic criteria for the determinacy of a germ /(x) are based on the so- 
called Jacobian ideal, J[/(x)]. This is the ideal generated by df/dxi, i.e. 
J[/(x)] = {gidf /dxi + . . . + gndf /dxn} for arbitrary germs 5i(x). By means 
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Name 


k 


r 


cat. germ 


unfolding 


fold 


A2 


3 


1 




aix 


cusp 


^3 


4 


2 


±x^ 


aix + a2X 


swallowtail 


Ai 


5 


3 




aix + a2X^ + 032;'^ 


butterfly 


^5 


6 


4 


±x'^ 


aix + a2X^ + a^x^ + a4X* 


wigwam 


Ae 


7 


5 


x' 


aix + a2a;'^ + asx'^ + a4x'* + 050;'' 



Table 1: Thorn's elementary catastrophes of corank 1. k: determinacy of germs; r 
codimension 



of the Jacobian ideal, we may also define the codimension of a germ /(x) as the 
dimension of the quotient space 

(3) cod(/) = dim(m„/J[/]) . 

We will go over the concepts of determinacy and codimension by discussing the 
Mather's criterion in Section [31 



3 Criterion for tricritical points 

Here we will assume that we are given a system whose phases are described by 4 
order parameters (namely x^y, z,w) and that the equilibrium phase is attained 
at a minimum of a free energy potential Vixjij^z^w) which is taken to be a 
smooth function. Moreover, we assume that V also depends on A; < 5 physical 
parameters, i.e., at the most 4 model parameters Xt plus the temperature T. 
For ease of notation, the physical parameters are not always explicitly indicated, 
but it is understood that V depends on them. As a consequence, also the 
position of the critical points will depend on Ai and T. We will follow [9] 
and define a tricritical point as a point in the phase diagram where a first 
order transition becomes second order. Our aim is to give, from a catastrophe- 
theoretic standpoint, a criterion that the physical parameters have to satisfy in 
order to identify a tricritical point in the phase diagrams of complex systems 
such as smectic or biaxial liquid crystals. 

In the simplified model describing the phase of a system in terms of a single order 
parameter ?/;, the free energy is usually given in terms of a Landau expansion as 

(4) V ^ + ia4V/ + 502^'' 

where 04,02 are the control parameters. It is well known [31 [S] that second 
order phase transition points in the phase diagram are given by 02 = 0, 04 > 
and the existence of a tricritical point is determined by a2 — = 0. This is 
easily understood by drawing the family of potentials ^ as 02 and 04 are varied 
(see Figure [T|) . The first condition corresponds to the appearance of two new 
minima placed symmetrically with respect to ?A = 0, while the phase V' = 
becomes unstable and the transition to the new phase is therefore continuous or 
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second-order. The second condition marks the transition to a higher degenerate 
potential V where the phase at -!/> = is still stable, therefore from there on the 

phase transition tn a. ih is pvnprtpA t.n hp first nrdpr 



X 



Figure 1: Potential @ as the control parameters are varied, (a) 02 > 0, 04 > ; (b) 
02 < 0, 04 > ; (c) a2 > 0, 04 < 0. 

We would like to extend the simple 1-variable tricriticality criterion a2 = 
04 = to more complex cases, preserving at the same time its valuable feature 
of being intuitively clear. Here, we have assumed the potential to be even, but 
this condition will be dropped in the following general treatment. A symmetric 
case will be considered in Section IH 

Of course, it is not possible to simply say "the fourth coefficient of the Taylor 
expansion at the transition is nought" as (a) the critical points arc not fixed 
but a slight variation of the control parameters implies a change in the equilib- 
rium order parameters and (b) there are more than one variable involved in the 
transition and it is not a-priori known along which direction the bifurcation is 
going to take place. In [531 [21] point (a) has been addressed by the introduction 
of a "master" order parameter which is different from zero only in the ordered 
phase. All the remaining order parameters are supposed to depend on it. We 
will not need to make this assumption here. 

At equilibrium, for given Xi and T, the order parameters solve W{x, y, z, w; A.;, T) = 
and we denote by Pq = (a:o(Ai, T), ?/o(Ai, T), zo(Ai, T), wo(Ai, T)) the critical 
point of our interest. To avoid unnecessary generality and to make our study 
close to the application in liquid crystal theory, we make the following assump- 
tions: 

1. V{x,y,z,w) has a corank 1 singularity at the transition. 

This means that the hessian matrix H is singular, i.e. det(]HI) — 0, but 
only one eigenvalue of H will be nought. Hence, we lose no generality in 
assuming that Vzz 7^ 

2. The third and fourth column of H are linearly dependent. 

Since from previous point det(H) = 0, we know that the columns of H 
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must be linearly dependent. Here, we want to express the condition that 
the bifurcation will involve only the z and w variables. After the transi- 
tion new minima will appear and if we check their location in the {z,w) 
plane, they will move along a direction determined by the eigenvector as- 
sociated with the null eigenvalue (see [4]). No bifurcation will involve x or 
y. Therefore, another way of stating the assumption is: the eigenvector 
associated with the null eigenvalue is of the form (0, 0, W3, V4). 

Now CT can be straightforwardly applied. Following assumption 1, the split- 
ting lemma allows us to study the functions with only one degenerate variable. 
Therefore, apart from a diffeomorphic change of variables that brings the "phys- 
ical variables" {x,y,z,w) into the "mathematical variables" (ai, 012, 0:3, "0) and 
the "physical parameters" (Ai,T) into the "mathematical parameters" a^; the 
potential V will locally be of the type: V{ai,a2, as, ip) = \a\ + -I- \a\ + 
f{ip]ai), where in the hst of elementary catastrophe (see Table[T]or the 
least degenerate (lowest codimension) example of a germ which is physically 
acceptable (bounded from below) and allows a tricritical point, is the butterfly 
A^: /{ijj) = ^tp^ + ja4^'^ + ^asip'^ + ^a2ij'^ + ai^ . In the case of even potentials, 
odd powers have to be ruled out for symmetry reasons. This is equivalent to 
the usual expression found in physics books [3]. 

By standard results in CT, all possible perturbations of the potential V in an 
open set around the critical point, i.e. all the possible topologically equivalent 
phase diagrams are given by the above expression when the at are varied. If, 
for ease of calculation, we assume that V has a critical point at -0 = 0, then 
ai — 0. It can now be shown that also in the multi-variable case, the phase 
tp — looses stability when 02 = and a second order transition happens at 
tp = when 02 = 03 = 0, a4 > 0. The presence of a term 03 ^ in general 
signals a first order transition to a phase with -0 7^ and therefore makes the 
whole study of the transitions much richer (and complex). Anyway, since we 
are here only interested in tricritical points, i.e. limit points where a transition 
pass from second order to first order, we do not need to exploit such complexity 
but can limit ourselves to the regions of the control space where second order 
transitions are possible. It is therefore recognised that the tricritical point in 
A'^ is given by what is konwn as the germ of the catastrophe, the point with 
maximum degeneracy (see [51 IH [HI [23] ) > given by 

(5) ai — a2 — — — 0. 

The condition ([5]) needs to be translated in terms of the physical parameters 
(Ai, T) if it has to be directly applied to a real world example, i.e., to a potential 
V{x,y, Zjw). This translation can be gained if we use the geometric language 
of CT to recast ([5]) in a intrinsic form, invariant over change of variables: (a) 
when the control parameters belong to the bifurcation set, then V{x, y, z, w) is 
not 2-determined and (b) when they belong to the tricritical set, V{x,y,z,w) 
is not 4-determined. In particular the necessary condition (b) for tricriticality 
is equivalent to (O and becomes also sufficient when it is not trivial, i.e., in the 
case of least codimension (butterfly A'^ catastrophe) considered above. 
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We thus arrive at the natural multi-variable extension of the criterion ^ on 
tricritcality: a necessary condition for the control parameters {Xi,T) to belong 
to the tricritical set is that the free- energy function V is not 4- determined at the 
critical point Pq (Ai , T) . 

The effectiveness of this approach lies in the powerful tools that are now at 
hand. In particular, Mather [51 Ej has given a completely algebraic procedure 
to calculate the determinacy of a function. 

Theorem 1 (Mather) Let V be a germ and let r > be an integer. If 

(6) m:+icmy[y]+m;+2, 

then V is r-determined. 

As usual, the sum and product of two ideals / and J are the two ideals 
defined by / + J := {i + j : « G /, j G J} and IJ := [ij : i e IJ £ J}. If I 
is generated by {ia} and J is generated by {jb}, then IJ is generated by the 
products iajb- Thus, we can define the powers 1^,1^, .. . of an ideal / taking 
the powers of its generators. The ideals mj'^ are indeed defined in this way as 
powers of m„. 

A condition stronger than ([6|), yet easier to check is 

(7) m;;-ic j[y] + m;, 

which follows form the properties of the ideals involved [5]. Then, if V is not 
4-determined the two conditions ^ and ([7]) with r ~ 2,3,4 have to be false. 
We will now translate the stronger expression ^ from the abstract form into 
an algorithmic form and apply it to our case. As usual in CT, in the abstract 
description we assume with no loss of generality that the Taylor expansions are 
performed in the origin. We will then express the results in terms of the poten- 
tial V evaluated in the critical point Pq (which is not necessarily the origin). 

Let us define Ai^-^ = j^To^ the r-jet linear space of the germs in m^. It 
contains all the linear combinations of the homogeneous monomials of degree 
r. To make our argument more concrete, we first concentrate on the case r = 3 
(the r = 2 case simply translates the bifurcation condition). Thus, we look for 
a necessary condition that non 3-determined functions must satisfy. The case 
r — 4:, which renders explicit the non 4-determinacy condition, is conceptually 
analogous and is obtained in exactly the same way. Taking the {i — l)-jct of ([7]) 
with r = 3 we have 

(8) MlcfUiV]} 

which of course needs to be false if V is not 3-determined. 
It is now easy to build a set of generators for {J[V^]}. It will be sufficient 
to give an example in the 1 variable case. Multi-variable case is conceptually 
analogous but simply involve cumbersome notations. Elements of J[V] are of 
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the type g{x)^, with g{x) a generic germ. We now suppose to expand g{x) in 
Taylor series up to the second order so that the remainder R{x) will be in mf^. 
Taking the 2-jet we have {gj £ K) 

(9) f {g{x) ^} = f{{go + 9iX + g2x' + R{x)) ^ } 

Therefore {J[V]} is generated by aU Uncar combinations of the polyno- 
mials of the type 9ij(x) = |pi(x)^^| where Pi(x) are all the monomials of 

degree < 2. Following [4j, we are now able to state Mather's theorem (in its 
strongest version ([7])) in algorithmic form. 

1. Calculate all the polynomials gij(x) — j^^^ |pi(x)^^|, where Pi(x) are 
all the multivariate monomials in n variables, of degree from to r — 1. 

2. If all the monomials of degree r — 1 can be obtained as linear combinations 
of the (7jj(x), then the potential V is r-determined. 

In our case {n — 4), when r = 3 there are 10 monomials of degree 2 and 20 
polynomials qij (x) whereas when r — A there are 20 monomials of degree 3 and 
60 polynomials qij (x). The number of the polynomials has been already reduced 
considering the fact that we take the jet in a critical point of V. However, 
further reductions are possible if we assume from the beginning that V is not 2- 
determined which implies that detH = (i.e. we want to study the transition). 
The number of zero terms is greatly enhanced when symmetry conditions are 
further assumed on V (see Section [5]). 

The monomials of degree from to 2 form a basis for both Ai^^ and {J[V]}j 
and thus the problem of 3-determinacy is reduced with respect to this basis to 
the solution of a set of linear systems Ac^ = b^, where the matrix A is built 
from the polynomials Qij , and 7 is an index that ranges only on the monomials 
of maximum degree. 

We assume here that we have chosen a degree term ordering within the mul- 
tivariate monomials of degree from to r — 1. The monomials of higher degree 
will be placed last with respect to this ordering. If there are Np monomials 
Pi and Nq polynomials gy , the matrix A and the vectors b^,c^ are defined as 
follows: 

• A e M(jVp,Arg) (K) is such that the entry in row s and column t is the coef- 
ficient of the s*^ monomial Ps (x) in the t*'* polynomial qij{^)- Therefore, 
along the columns of A we can read all the coefficients of the polynomial 
qij (x) , the s*'^-row being associated with the monomial ps (x) , 

• are column vectors in M(jVp.i)(I^) with all first zero elements and only 
a 1 in 7 position to represent one of the higher degree monomials. The 
index 7 will range in order to solve a linear system for each of such a 
monomial (i.e., x^y^ z^w*' such that i,j,h,k > and i + j + h + k = r — 1). 
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• c-y arc the vectors in M(jv .i)(R) eontaining the unknown coefficients that 
give the Unear dependence of the monomial Pj in terms of the polynomials 
qij- 

Here, M(/j is the space of real matrices with h rows and k columns. 

A necessary condition for non 3-determinacy is therefore that not all such 
linear systems can be solved. Since it is well known from linear algebra that 

e 7^(A) ^ b,,^ e {Af{A'^)) \ this means that our criterion on tricriticality 
requires that there exists a bo G A/'(A-^) such that bg • b-^ 7^ for at least one 
b-^. Here, 7^(A) is the range of A and J\f{A) is the null space of A. So we are 
now induced to solve the single system A-^bo = and see what condition on 
the entries of A can be imposed in order to get a non-trivial solution. In such a 
way we are simply trying to express one of the last rows of A (representing the 
monomials p-y) as a linear combination of all the other rows. This problem can 
be solved easily, only the dimensions of the matrix A suggests that it is wise to 
perform calculations with the aid of a symbolic computing software. The non 
4-determinacy is given by retracing the same algorithm with r = 4. The num- 
ber of polynomials involved is obviously raised and therefore the computations 
increase in complexity, but still a condition can be found. 

Let us define the matrix M as the 3 by 3 upper-left submatrix of H. It 
comprises the entries of H that are derivatives of V with respect to the variables 
X, y and z only. Let M = det(M) and Mij the (i,j)-minor of M, i.e., the 
determinant of the 2 by 2 matrix formed by removing from M its i*^ row and 
j*^ column. 

After some lengthy but easy algebra, we finally arrive at the conditions that a 
non 4-determined potential must satisfy. 

1. V is not 2-determined (and therefore a bifurcation occurs) if dct(]HI) = 0. 
By virtue of the assumptions on V it can be shown that the above equation 
is equivalent to the requirement that det(]Hl2) = 0, where H2 is the 2 by 2 
lower right submatrix of H. Therefore V is not 2-determined if 

(10) D2y = 0. 

2. V is not 3-determined if it is not 2-detcrmined and 

(11) D^V = 0. 

3. V is not 4-determined if it is not 3-determined and 

3 

(12) M - 3 ^ (-1)*+^ Mij KiKj = . 
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where 






(13) 


a 




(14) 


vP-v 


= cP'Vzz - "^olVzw + y^w , 


(15) 






(16) 


dV 


— V^2:2:2: ~ 4:0; \^zzzw H" So Vzzww ^ttVg^^^ -|- 'Vwwww 


(17) 




O VxZZ ' ^OV^^-LU H~ ^^lUlO ) 


(18) 


K2 


— CI Vyzz 2o;V^_2t(; H~ yyww ■> 


(19) 


K3 


— O Vzzz ^CkVg^^ -|- Vzww- 



A subscript denotes differentiation with respect to the indicated variables and 
all the expressions are evaluated at the critical point Po{Xi, T). We recall that 
according to our assumption 1 above, it is Vzz ^ 0. 

In conclusion, (fTUl) - P^ are three generally independent equations in the 
control parameters (variables are evaluated at the critical point). They together 
define the set of control parameters at which a tricritical point occurs in the 
phase diagram. More precisely, if we have k control parameters they define a 
fc-3 dimensional locus of tricritical points. 



4 Reduction to simpler cases 



There are important simplifications to the general result that suit many physical 
applications. We present them here, adding more constraints to the potential 
V as we go on in the Section. 

First, we assume that V satisfy an even symmetry with respect to a simultaneous 
change of sign (z, w) ^ {—z^ —w) 

(20) V{x,y,-z,-w) = V{x,y,z,w). 

This of course implies that at the critical point 

QP+q+r+sy 



(21) 



dxPdyidz^dw'' 



= 0, 



if' 



s is odd. 



The equilibria with z = w = are then the natural candidates for phases 
from which a second order transition could develop. Therefore, we assume the 
critical point under study is Pq = {x(){Xi,T)^yo{Xi,T)^0,0) which, as before, 
generally depends on the control parameters (Ai,T). Moreover, the hessian 
matrix will be block diagonal: 





V v^z 



\ 



Vz. 
^ww j 
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As in the general case, the bifurcation will involve only the z and w variables. 
Explicitly, in terms of the hessian matrix, we have that Hi is positive definite 
(det(]HIi) > and V^x > 0) and H2 has a null eigenvalue at the transition (and 
only one since we are assuming a corank 1 singularity). 

Furthermore, this symmetry directly implies that V is not 3-determined and 
PT|) is automatically satisfied. This condition is found in many physical 
applications where a second order transition occurs. Therefore, the reduced 
equations are worthy of being stated explicitly 

(22) K,K,^ - V,l = detH2 = , 

(23) {VxxVyy - Vly) D^F - 3 {k\ Vyy ~ Vxy + kI Vxx) = . 

The choice of 4 order parameters has been mainly motivated by the desire to 
apply our criterion to existing Landau free energy expansions in liquid crystals 
theory, in particular to smectic or biaxial liquid crystals. However, many simpli- 
fications or reduced Landau expansions can be given where the phase transitions 
are well described by fewer order parameters. Of course, ([22]) and ((23|) ought to 
reduce to the equations already existing in literature when we further simplify 
the free energy V . 

For instance, in the case of two order parameters, we can declare 2 out of 
4 order parameters to be "dummy variables" . This means that a minimization 
with respect to this variables would always yield a constant value (which can 
be assumed to be zero) and no bifurcation in these variables can occur. 

The free energy will therefore be 

(24) V{x, y, z, w) ^ iy2 ^ ^ ^(^^ ^) 

where, in agreement with our previous assumption on V, g{x, —w) — g{x, w) 
and at the critical point > 0. We remark here that g{x,w) depends also 
on the control parameters, even if not explicitly indicated. Equations (j22p and 
((231) yield 

(25) g^iw = 0, 
(26) 

where again the derivatives have to be evaluated at the critical point. The same 
result is to be found in [TH [13] (with x = S,y = S',z = T and w = T') where 
it is applied to biaxial liquid crystals in the fairly general yet simplified case of 
a null model parameter. Indeed an extension of that result has been one of the 
purposes of the present work and in fact (|22p and (|23p can in principle give the 
locus of tricritical points in the phase diagram of biaxial liquid crystals in their 
general treatment employing 4 order parameters. Unluckily, it is not at the 
present time possible to pursue this calculations any further because a reliable 
Landau expansion of the free energy, using 4 order parameters, is not to our best 
knowledge available yet. We are anyway aware that some work is being carried 
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on by few research groups [TU] and therefore we intend to perform a study of 
the tricritical locus in biaxial liquid crystals in a subsequent paper when these 
results are published. On the other hand, we could in principle apply ([^ and 
(j23p to a non polynomial expression, such as the integral representation of the 
free energy given by a mean field model [HI [TU [131 [HI [15] . Yet, this could be 
done only at the cost of great computational efforts and employing numerical 
approximations . 

As a last example we consider a system which is described just by one order 
parameter w. The free energy can be given the form 

(27) Vix, y, z, w) = Ix" + iy^ ^ i ^2 ^ j(^) 
where /(— w) = /(w). A substitution in (|^^ and (|23p yields 

(28) fww — fwwww — , 

in agreement with the standard condition ([5]). 

Of course, when necessary one could gather other possible reduced equations 
according to various combinations of dummy and non-dummy variables. The 
examples given here have the immediate advantage of reproducing the structure 
of the Landau expansions usually found in the physical literature. 



5 Application to smectic liquid crystals 

In this Section we apply our result to the case of smectic-C liquid crystals, 
where a tricritical point along the Sm-A/Sm-C transition line has been reported 
[H [m [TSl [HI HD]. We refer to [HI HD] and bibliography therein for the physi- 
cal details and to [50] for the expression of the Landau expansion we are going 
to use. Nematic liquid crystals [TJ [2] are formed by elongated molecules with a 
cylindrical symmetry. When the thermal motion is not dominant (at sufficiently 
low temperature), they tend to align themselves along a common direction, usu- 
ally identified by a unit vector n, called the director and form a nematic (N) 
phase. The degree of orientation is measured by an order parameter s (where 
s = and 5 = 1 correspond to no order and to perfect order, respectively). 
The smectic phase is established when an additional positional order is super- 
imposed to the orientational order: the centre of mass of the molecules are 
organized in layers. The layers are described by a periodic modulation of am- 
plitude p in the mass density of molecules in the liquid. An additional order 
parameter q gives the inverse of the distance of two subsequent layers: d = 2'K/q. 
When the director n is normal to the layers, the phase is called smectic-A (Sm- 
A). The onset of a smectic-C (Sm-C) phase is described by a non-zero angle 9 
between the director and the normal to the layers. We will use w — sin 9 as the 
order parameter for the Sm-C phase. 

It is known that some liquid crystal compounds show a first order Sm-A/Sm-C 
transition (for instance B7) while other undergo a second order transition (for 
instance 10O4). We follow [21] and study a binary mixture of two such com- 
pounds. We imagine to increase the concentration of one of the compounds 



A Catastrophe-Theoretic Approach to Tricritical Points 



13 



and alter the nature of the transition. It is then reached a concentration value 
where the phase diagram shows a tricritical point. In a homogeneous state, the 
polynomial expression for the free energy density of the binary mixture is given 
by expression (2.5) of [20] : 

(29) V = ias^ - i6s3 + i^^4 ^ i ^^2 ^ i^^4 ^ ig^ ^2 ^ 1 ^^2^2 

Here we depart slightly from the notation found in [20] not to generate confusion 
with the symbols we have already adopted in previous Sections with different 
meaning. The notation in [20] is fully restored when we put p — ipOi a — a and 

The order parameters are: s, p, q and w. They identify the various phases 
whenever assume a non-zero value. More precisely: isotropic phase (I) is s = 
p ^ q ~ w = 0; nematic phase (N) is s ^ 0, and p — q — w = 0; smectic-A 
phase (Sm-A) issy^O, pj^Oqy^O and w — 0; smectic-C (Sm-C) phase is 
s^O, p^O, q^O and w ^0. At a first sight one might expect the Sm- 
A/Sm-C transition to be only second order because the expansion is even with 
respect to the order parameter w and has been truncated at a point where w 
appears only to the fourth power. It must be noticed, however, that a first 
order transition is possible due to the coupling between w and the other order 
parameters. One way to see this, is to think w as a master order parameter and 
use the equilibrium equation Vs — 0, Vp — and = to express s, p and q as 
functions of w. It is then apparent that the various terms involving a coupling 
between s, p, q and w will bring higher terms in w into consideration, therefore 
giving an intuitive justification for the presence of a tricritical point. 

Following the ideas of McMillan and de Gennes [24] about the character 
of the N/Sm-A phase transitions, something can be argued about the transi- 
tion temperatures at the tricritical point. It is well known (as often said, due 
to presence of a cubic term) that s undergoes a first order transition I/N at 
the transition temperature T/v/. If Tag is the temperature at which the Sm- 
A/Sm-C transition takes place, close to the nematic phase (Tac/Tni ~ 1) the 
fluctuations of orientations are large, and as a result the character of the transi- 
tion is determined by the nematic order parameter s, i.e., the phase transition is 
first order. For low transition temperature {Tac/Tni ~ 0), the nematic order is 
nearly saturated and thus the character of the transition is driven by w and will 
be second order. Therefore, the tricritical point along the Sm-A/Sm-C transition 
line is expected to be in a region of the control space where Tag < Tma < T^i 
but Tag ~ T^a ~ T^i- 

Suppose now that we want to study the nature of the transitions from a Sm- 
A to a Sm-C phase. Minimization of ((^U) with respect to s, p, q and w yields 
the equilibrium values in terms of the control parameters. We do not explicitly 
write the expressions here, simply write the equilibrium order parameters in the 
Sm-A phase as s = s+ > 0, p = /?+ > 0, g = (7+ > and w ^ Q. With the 
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identification x = s, y — p, z = q and w = w, we can use the piI)) - ([T^ to find 
the conditions that the control parameters must satisfy in order to go from a 
second order to a first order Sm-A/Sm-C phase transition. 
The bifurcation condition (fTO|) yields 

2 

+ ' 



(30) e = -4/i s+ q 



while ([TT|) is trivially satisfied. 

Equation (|30p is not sufficient to guarantee that a transition takes place, as the 
transition might be first order. The limit point of second to first order transition 
is yielded by p2|) which, after a little algebra, reads 



(31) 

_ v° (V° ? - (V° ?V° + 2V° V° V° 

ss\ pq / \ sq / pp sp sq pq 



4 



where 



(32) V°, =a- 2bs+ + 3csl + jpl + ^hp\qX + rj^ + 2ujs+^ , 

(33) T/° ^{S + 2js+)p+ , 

(34) V^^ =8hs+plql, 

(35) V°p =a + X + Ss+ + 74 + 3/3pl + diq\ + ^2?+ - 4/i44 , 

(36) -(2di+4d2g^)p+, 

(37) = (di + 6^29^ + 16/^44)4 , 



are the entries of the hessian matrix evaluated at the critical point (s_|_, p+, (7+, 0) 
and (|30l) has been used. 

Equation (j31|) can now be solved with respect to any of the control parameters 
(paying attention that also the equilibrium values s+, p+ and q+ depend on 
the control parameters). Usually in Landau expansions it is assumed that only 
the coefhcients of second order terms depend on temperature, linearly. This is 
usually the first step towards a physical interpretation of (l3T|) . For instance, we 
can take a — ao{T — Tj*) and a = do(T' — T2*), where T is the temperature, Tj* 
is the nematic supercooling temperature and T2* is the virtual (second order) 
transition temperature into the smectic phase. At Tj* the solution s = s+, p = 
looses stability. The tricritical temperature is then given by solving ||3T|) with 
respect to the temperature T. It is remarkable that the same results found in 
[20] after quite lengthy and tedious calculations, are given directly by ([3T|) . In 
particular, ((3T|) is trivially satisfied when we substitute all the expressions that 
lead to the tricritical temperature (2.23) and the tricritical concentration (2.24) 
in [20j . However, physical significance must be given to the control parameters, 
for example by means of a comparison with experimental results, before physical 
conclusions could be drawn. This is not in general an easy task and it is outside 
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the scopes of the present paper. 

We must note that in general ([3T|) does not guarantee that a tricritical point is 
found, for the simple fact that the critical point we are exploring might not be 
the absolute minimum of the free energy. Our analysis is local. The event that 
the non 4-determinacy happens in a local minimum which is not an absolute 
minimum must be studied by means of other global techniques and a more 
comprehensive description of the whole phase diagram is needed. 
Finally, we would like here to underline that our theory can in principle be 
applied to more complex potentials, such as a free energy functional derived 
by a mean-field model (see for instance [12l [l3j for example of applications of 
(P5l) . (|26p to Sm-A and biaxial liquid crystals). 

6 Conclusions 

General criteria usually found in literature to locate tricritical points in the 
phase diagram loose their elementary interpretation when applied to complex 
systems such as smectic or biaxial liquid crystals. Moreover, they might not be 
applicable to these materials when the full generality on the number of physical 
parameters is allowed. 

We have shown that rephrasing the simple and intuitive criterion employed for 
single order parameter systems using the invariant language of CT, allows an 
immediate extension to the multi order parameter cases. Therefore it is possible 
and often easy to explicitly write the equations that identify the tricritical points 
for instance in the rich and complex scenario of smectic. This extension was 
in fact the main aim of our work. We believe that also the interesting scenario 
of biaxial liquid crystals, within the general description which comprises four 
order parameters, can be addressed with the results described in the present 
paper. Our criterion can thus help in giving an analytic answer to the final 
question posed in [12]. There the authors raise the issue whether the predicted 
tricritical point at the direct transition between the isotropic and the biaxial 
phases persists and could possibly extend when all physical parameters (and 
consequently all four order parameters) are taken into account. Indeed, we 
aim at undertaking a more comprehensive study of this problem, relying on a 
Landau expansion of the free energy which seems not available at the moment 
(CHECK), in a subsequent paper. 
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